Probing the Strong Boundary Shape Dependence of the Casimir Force 



O 

o 

(N 

> 
O 

00 



o 

B 



X3 

o 

> 
oo 

O 

O 



I 

o 
o 



X 



Thorsten Emig^, Andreas Hanke^, Ramin Golestanian^^"', Mehran Kardar^^'^ 
^Physics Department, Massachusetts Institute of Technology, Cambridge, MA 02139 

"^Institute for Advanced Studies in Basic Sciences, Zanjan 4.5195-159, Iran 
^Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106 
"^Institute for Studies in Theoretical Physics and Mathematics, P.O. Box 19395-5531, Tehran, Iran 

(February 6, 2008) 

We study the geometry dependence of the Casimir energy for deformed metal plates by a path 
integral quantization of the electromagnetic field. For the first time, we give a complete analytical 
result for the deformation induced change in Casimir energy 5£ in an experimentally testable, 
nontrivial geometry, consisting of a flat and a corrugated plate. Our results show an interesting 
crossover for 5£ as a function of the ratio of the mean plate distance H, to the corrugation length 
A: For \ <^ H we find a slower decay ~ i?"**, compared to the behavior predicted by the 

commonly used pairwise summation of van der Waals forces, which is valid only for H. 

PACS numbers: 03.70.+k, ll.lO.-z, 42.50.Ct, 12.20.-m 



In 1948, Casimir showed that two parallel conducting 
plates, separated by a distance H, attract each other with 
a force F, proportional to the surface area A, and given 

by i 
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This remarkable prediction of quantum electrodynamics 
can be understood as resulting from the modification of 
the zero point vacuum fluctuations of the electromagnetic 
field by the presence of boundaries. Since this discovery, 
the fundamental nature of the Casimir effect and its im- 
plications, e.g., on surface forces ||], particle physics 
and cosmology |4|, has motivated extensive theoretical 
work (see Refs. p-M for reviews). On the experimental 
front, the initial attempt at observing the Casimir force, 
by Sparnaay in 1958, was not conclusive due to large ex- 
perimental uncertainty jl^. Only recently, there have 
been a number of precision measurements of the Casimir 
force, using a torsion pendulum \ ^a\ , an atomic force mi- 
croscope [ p^ , and a micromachined torsional device , 
which confirm the theory to a few per cent accuracy. The 
latter experiment also demonstrates the possibility for 
novel actuation schemes in microelectromechanical sys- 
tems based on the Casimir force [pif . 

In the more general context of the Lifshitz theory for 
dielectric bodies |l|], Eq. @) appears in the limit of per- 
fectly conducting plates, for which the dielectric constant 
e is infinite. For finite e = s{uj), this power law for the 
force is recovered for large distances H ^ c/ujq, where 
luq is the smallest resonance (absorption) frequency of the 
dielectric (usually c/ujo « 10-100 nm). In this, so-called 
retarded, limit, the force is universal in the sense that 
it only depends on the electrostatic dielectric constant 
Eq — e(0), and can be obtained, e.g., by dispersion re- 
lation techniques |^6|. The opposite limit of H ^ c/luq 
gives the unretarded van der Waals force F/A ~ H"^, 
which can also be obtained by summing the (attrac- 



tive) intermolecular interactions due to induced molec- 
ular dipole moments. Even though obtainable from the 
same microscopic theory, the Casimir and van der Waals 
forces are quite different. In particular, the interpreta- 
tion of the Casimir force in terms of changes in zero point 
vacuum electromagnetic energy suggests it to be a strong 
fimction of geometry |l^,|l^; probing the global shape 
of the boundary that confines the vacuum fluctuations. 
Indeed, whereas the van der Waals force between elec- 
trically polarizable particles is always attractive, even 
the sign of the Casimir force is geometry dependent, 
and can be repulsive, e.g., for a thin spherical or cubic 
shell [^,@jl^] . (Repulsive Casimir forces are expected also 
when magnetic as well as electric properties are included 
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FIG. 1. Geometry used for calculating the Casimir energy 
of a flat plate and a corrugated plate at mean distance H. 

It is highly desirable to demonstrate the strong shape 
dependence of the Casimir interaction in a set-up that 
clearly demonstrates its distinction from the usual pair- 
wise additive interactions pof . Since measurement of 
the repulsive Casimir interaction for a conducting sphere 
is experimentally difficult, the most promising route is 
to via modifications of the parallel plate geometry. In 
searching for nontrivial boundary dependences, Roy and 
Mohideen ||2^ examined the force between a sphere, and 
a sinusoidally corrugated plate with amplitude a « 60 nm 
and wavelength A ~ 1.1 /im. Over the range of separa- 
tions H fa 0.1 — 0.9 /Ltm, the observed force showed clear 
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deviations from the dependence expected on the basis 
of decomposing the Casimir force to a sum of pairwise 
contributions (in effect, an average over the variations 
in separations). This experimental result motivated our 
calculation of the exact Casimir force in the geometry 
depicted in Fig.|l|, without the assumption of pairwise 
additivity. Our analytic results [see Eqs. (4) - (6) and 
Fig. H hold to second order in a, and show that for fixed 
H the corrections due to corrugation strongly depend on 
A. In fact, for H/\ ^ 1 the correction is by a factor 
of H/X larger than in the opposite limit of H/X <C 1 
where the assumption of pairwise additivity is asymptot- 
ically correct. However, the experiments of Ref. are 
performed in the range of H/X k, 0.1 — 0.8 where the 
corrections to pairwise additivity are in fact not signifi- 
cant enough to account for the observed deviations. This 
bolsters the conclusion in Ref. that these deviations 
are in fact due to a lateral force that tends to prefer- 
entially position the spherical AFM tip on top of local 
maxima of the modulated surface (leading to a smaller 
separation and stronger force). We thus propose that 
the shape dependence of the Casimir force can in fact be 
probed in this set-up by going to modulations of shorter 
wavelength; a hard but achievable goal. 

The use of a spherical tip, of large radius R, in experi- 
ments [pT|-p^ causes some differences from the flat plate 
geometry used in our calculations. First, the position- 
ing of the tip relative to the modulations is important 
when H and A are comparable, but becomes insignifi- 
cant in the proposed limit of A <C i?, R- Secondly, as 
long sls R ^ H,X the curvature of the tip does not lead 
to nontrivial corrections, and the force can be related to 
the energy per surface area £ in Eq. (5) by the proximity 
force rule F = 2ttR£ . These formulea thus provide a 
specific recipe for evaluating the nontrivial shape depen- 
dences of the Casimir force in the experimental set-up. 

Before turning to the geometry of Fig. |^, consider the 
more general case of two perfectly conducting plates of 
mean separation H, which are infinitely extended along 
the x-y plane. Now assume that one of the plates is de- 
formed in a way that is translationally invariant along the 
y axis and has no overhangs. Its profile can then be de- 
scribed by the height function h{x), with J dxh{x) = 0. 
For example, to describe the geometry in the experiment 
of Ref. Il^, we choose h{x) = acos(27ra:/A) as in Fig. |^. 
The Casimir energy associated with general h{x) at zero 
temperature corresponds to the difference of the ground 
state energies of the quantized electromagnetic (EM) field 
with and without plates, respectively. To obtain this en- 
ergy, we employ the path integral quantization method, 
which can be applied to the EM gauge field by intro- 
ducing a suitable gauge fixing procedure j24|. However, 
in the present translationally invariant geometry we can 
develop a simpler quantization scheme, by a similar rea- 
soning as used in the context of waveguides with constant 



cross-sectional shape (here along the y axis) . For an 
arbitrary EM field between the plates, transverse com- 
ponents of the E and B fields are completely determined 
by their axial components Ey and By. Therefore, any 
EM field can be described by a superposition of two in- 
dependent scalar fields $tm = Ey (transverse magnetic 
waves. By = 0) and $te = By (transverse electric waves, 
Ey — 0). The scalar fields both fulfill the usual wave 
equation, but differ in their boundary conditions on the 
plates S, as <I>tm|s = 0, while 9n<&TE|s = 0, where 9„ 
denotes the normal derivative. 

Both scalar fields can now be quantized by considering 
the Euclidean action S'[(i>] = 5 / d'^X (V<i>)^, correspond- 
ing to the wave equation after a Wick rotation to the 
imaginary time X'^ = ict. In the 4D Euclidean space, 
the plates are parametrized by Xi{r) = [r,h{x)] and 
^2(r) = [r, ^f], with r — {ict,x,y). We implement the 
boundary conditions on S using delta functions |p0| , p6| , 
leading to the partition function 

2 

h'^il C[^X,)] eM-Sm/h), (2) 

with the boundary condition enforcing functionals 
CMX,)] = YljmXjir)) for $tm, and C[$(X,)] - 
rir ^{'^ri^{Xj{T)) for <I>TE, and the partition function Zq 
of the space without plates. The Casimir energy per sur- 
face area A is then given hy £ = —hc\nZ/AL where L 
is the overall Euclidean length in time direction. Imple- 
menting the delta functions by integrals over auxiliary 
fields and integrating out $, we obtain an h{x) depen- 
dent kernel M for the Gaussian action of the auxiliary 
fields. Expanding Z — (detAf)"^/^ to second order in 
h{x), we get for the deformation dependent part of Z, 

\nZh^^ IJ^ K{r-r')h{x)h{x'). (3) 

The kernel is the sum of contributions from the two 
wave types, i.e., K{r) = ifxividrl) + isrTE(koUr|| I), with 
r = (ro,r||), and has been calculated explicitly. 

For the specific deformation of the plates correspond- 
ing to harmonic corrugation of amplitude a and wave- 
length A defined above, the calculation of In Z^ reduces to 
Fourier transforming the kernel K{r). The correspond- 
ing integrals can be performed for A > by closing the 
integration contour via a semi-circle at infinity in the up- 
per half of the complex plane. The resulting sum of an 
infinite series of residues can be expressed in terms of the 
polylogarithm function Li„(z) = X^l^i ^"^/f^", leading to 

^ = -W^{y20 + 1p [^™ (f) + ^™ (f)]} ■ 

(4) 

The contributions from TM and TE modes are: 
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with u = exp(— 47rs). Figure ^ displays separately the 
contributions from Gtm and Gte to the corrugation 
induced correction S£ to the Casimir energy. While 
Gtm{H/X) is a monotonically increasing function of 
H/X, Gte{H/X) displays a minimum for i// a w 0.3. The 
net Casimir energy £ is shown in Fig. ^ for two represen- 
tative values of a/A, including the parameters used in 
the experiment of Ref. |2^. Note that the corrugation 
induced correction leads to a larger energy £, and hence 
the corresponding force F — 2ttR£ is enhanced. 
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FIG. 2. Rescaled correction 5£ to the Casimir energy due 
to the corrugation as given by the terms in square brackets 
of Eq. (upper curve) . The lower curves show the separate 
contributions from TM and TE modes. The rescaling of S£ 
is chosen such that the result from a pairwise summation of 
van der Waals forces is a constant (dashed lines). 

Examining the the limiting behaviors of Eq. is in- 
structive. In the limit X ^ H, the functions Gtm and 
Gte approach constant values, and the Casimir energy 
takes the A-independcnt form 

he TT^ f 

IP 720 V + ' IP )■ 

Note that only in this case both wave types provide the 
same contribution to the total energy, see Fig. In the 
opposite limit of A -C -ff , as suggested by the first terms 
in Eqs. (||), (^, both Gtm and Gte grow hnearly in H/X. 
Therefore, in this limit the correction to the Casimir en- 
ergy decays slower, according to a new power law in H, 

he TT^ / \ 
£ = H-27r , (8) 



with an amplitude proportional to 1/A. Analyzing the 
correction S£ in the limit a, X <^ H for arbitrary val- 
ues of a/ A, we find that the factor multiplying a/H in 
Eq. (^) saturates for A <C a at a number of order unity. 
This result can be justified by noting that the most rele- 
vant contributions to the force come from modes of wave- 
length of order H. The corrugation also affects modes of 
wavelength of order A, but these modes contribute to 
the single plate energy only. Thus, in the extreme limit 
A -C a, one has a clear separation of the length scales 
H and A, and the modes "see" flat plates at an effective 
separation H — a, leading to a correction of the order 
a/H after expansion in a. 
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FIG. 3. Rescaled Casimir energy as given by Eq. (M) for 
two fixed values of a/A. The rescaling is chosen such that 
the Casimir energy of two flat plates becomes one (horizontal 
line). The lower curve with a/A — 0.05 corresponds to the 
parameters used in the experiment of Ref. IT), where H/a 
varies between approximately 3 and 17. 

The above behavior of the correction S£ for small and 
large H/X clarifies the limits of validity of previous re- 
sults in the literature. The upper dashed line in Fig. ^ 
corresponds to a widely used approach P,p7t in which 
the interaction is obtained from a pairwise summation of 
'van der Waals type' two body forces. It is evident that 
this approximation is accurate only for H/X — > 0, which 
in this limit is equivalent to the Derjaguin method to any 
order in the amplitude a p^ . Already for H/X of order 
unity, the additive van der Waals type approximation 
breaks down. The opposite limit, H/X — > oo, corrobo- 
rates the result reported in Ref. , which is larger than 
the former by a factor of H/X ^ 1. However, in experi- 
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ments with lateral distortions A of the order of iJ, none 
of the above limiting cases is realized, which makes the 
present, more complete analysis necessary. 

Moreover, for purposes of experimental comparison, 
corrections due to finite conductivity of the plates, sur- 
face roughness, and finite temperature should be taken 
into account. These corrections introduce additional 
length scales into the problem, which are in turn the 
plasma wavelength Xp of the plates (e.g., Xp ~ 100 nm 
for aluminium the transverse correlation length 

^ of the roughness (usually ^ « 300 nm |13|), and the 
thermal wavelength Xt = hc/ksT (« 1 /im at 300° K). 
The plasma and thermal wavelengths provide lower and 
upper bounds for H, respectively, such that our results 
for perfectly conducting plates at zero temperature are 
vahd for A, iJ » Ap, and < At 

Finally, we note that in the set-up of Fig. || nontriv- 
ial shape dependencies appear as corrections to a larger 
Casimir force. For the purpose of experimental tests, it 
is much more desirable to devise set-ups which directly 
probe differences, without the need for subtracting a 
larger baseline force. For example, in an atomic force ex- 
periment, simultaneous scanning of a flat and corrugated 
substrate would be desirable; while in the torsion pendu- 
lum experiment, one can imagine suspending a spherical 
lens equidistantly from two plates, one of which is cor- 
rugated. Another potential experiment along these lines 
is to measure the lateral force between two plates with 
sinusoidal corrugations of the same wavelength A, which 
are shifted relative to each other by a distance S. We 
find a lateral force F|| = fi-c^ff^ sm{2TT S / X) g{H / X) . This 
force tends to change the position of the plates such that 
a maximum is opposite to a minimum, corresponding to 
6 — A/2. The universal function g{H/X) tends to a finite 
value for H/X — > 0, but vanishes exponentially for large 
H/X. Similar results hold for the alternating component 
of the standard Casimir force as a function of the phase 
shift between the plates. 

We thank J.-P. Bouchaud, S. Dietrich, S. G. Johnson, 
M. L. Povinelli, and S. Scheidl for useful discussions. This 
work was supported by the Deutsche Forschungsgemein- 
schaft under grants No. EM70/1-3 (T.E.), HA3030/1-2 
(A.H.), and the National Science Foundation through 
grants No. DMR-01-18213, and PHY99-07949 (M.K.). 



[1] H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 
(1948). 

[2] J. N. Israelachvili, Intermolecular and Surface Forces 

(Academic, London, 1992). 
[3] K. A. Milton, Phys. Rev. D 22, 1441; 22, 1444 (1980). 
[4] I. Brevik and H. Kolbenstvedt, Nuovo Cimento B 82, 71 

(1984). 

[5] G. Plunien, B. Miiller, and W. Greiner, Phys. Rep. 134, 
87 (1986). 



[6] V. M. Mostepanenko and N. N. Trunov, The Casimir 
Effect and its Applications (Clarendon, Oxford, 1997). 

[7] E. Elizalde and A. Romeo, Am. J. Phys. 59, 711 (1991). 

[8] P. W. Milonni, The Quantum Vacuum (Academic, San 
Diego, 1994). 

[9] M. Kardar and R. Golestanian, Rev. Mod. Phys. 71, 1233 
(1999). 

[10] M. J. Sparnaay, Physica (Utrecht) 24, 751 (1958). 

[11] S. K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997); 81, 

5475(E) (1998). 
[12] U. Mohideen and A. Roy, Phys. Rev. Lett. 81, 4549 

(1998) . 

[13] H. B. Chan, V. A. Aksyuk, R. N. Kleiman, D. J. Bishop, 
and F. Capasso, Science 291, 1941 (2001). 

[14] F. M. Serry, D. Walliser, and G. J. Maclay, J. Microelec- 
tromech. Syst. 4, 193 (1995). 

[15] E. M. Lifshitz, Sov. Phys. JETP 2, 73 (1956); L E. 
Dzyaloshinskii, E. M. Lifshitz, and L. P. Pitaevskii, Adv. 
Phys. 10, 165 (1961). 

[16] G. Feinberg and J. Sucher, Phys. Rev. A 2, 2395 (1970). 

[17] R. Balian and B. Duplantier, Ann. Phys. (New York) 
104, 300 (1977); 112, 165 (1978). 

[18] However, one should note that by virtue of Schwinger's 
source theory [J. Schwinger, L. L. DeRaad, Jr., and K. 
A. Milton, Ann. Phys. (New York) 115, 1 (1978)], the 
Casimir force can be obtained without referring to vac- 
uum fluctuations; see also Ref. 

[19] T. Boyer, Phys. Rev. A 9, 2078 (1974); V. Hushwater, 
Am. J. Phys. 65, 381 (1997). 

[20] R. Golestanian and M. Kardar, Phys. Rev. Lett. 78, 3421 
(1997); Phys. Rev. A 58, 1713 (1998). 

[21] A. Roy and U. Mohideen, Phys. Rev. Lett. 82, 4380 

(1999) . 

[22] G. L. Klimchitskaya, S. L Zanette, and A. O. Caride, 

Phys. Rev. A 63, 14101 (2000). 
[23] B. Derjaguin, Kolloid Z. 69, 155 (1934). 
[24] M. E. Peskin, and D. V. Schroeder, An Introduction to 

Quantum Field Theory (Addison- Wesley, Reading, 1995) . 
[25] J. D. Jackson, Classical Electrodynamics (Wiley, New 

York, 1999). 

[26] Following the introduction of the path integral method, 
further apphcations have been considered in O. Kenneth 
and S. Nussinov, Phys. Rev. D 63, 121701 (2001); for a 
single corrugated plate, see also A. Hanke and M. Kardar, 
Phys. Rev. Lett. 86, 4596 (2001). 

[27] M. Bordag, G. L. KUmchitskaya, and V. M. Mostepa- 
nenko, Mod. Phys. Lett. A 9, 2515 (1994); Int. J. Mod. 
Phys. A 10, 2661 (1995). 

[28] S. K. Karepanov, M. Y. Novikov, and A. S. Sorin, Nuovo 
Cimento B 100, 411 (1987). 

[29] The importance of stochastic surface roughness can also 
be deduced from our calculations. The relative correc- 
tions 5£/£ to the Casimir energy due to roughness of 
amplitude a and transverse correlation length ^ should be 
of the form / H'^ for H, and a'^/{^H) for ^ < //. 
While the latter behavior is in accordance with Ref. [^ , 
we note that the experimental case corresponds to neither 
extreme, making a more complete analysis necessary. 

[30] M. Y. Novikov, A. S. Sorin, and V. Y. Chernyak, Theor. 
Math. Phys. 82, 124 (1990); 82, 252 (1990); 91, 658 
(1992); 92, 773 (1992). 



4 



